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Abstract. Preliminary calculations using the Covariant Spectator Theory (CST) em¬ 
ployed a scalar linear confining interaction and an additional constant vector potential to 
compute the mesonic mass spectra. In this work we generalize the confining interaction 
to include more general structures, in particular a vector and also a pseudoscalar part, as 
suggested by a recent study [I|]. A one-gluon-exchange kernel is also implemented to 
describe the short-range part of the interaction. We solve the simplest CST approxima¬ 
tion to the complete Bethe-Salpeter equation, the one-channel spectator equation, using 
a numerical technique that eliminates all singularities from the kernel. The parameters of 
the model are determined through a fit to the experimental pseudoscalar meson spectra, 
with a good agreement for both quarkonia and heavy-light states. 


1 Introduction 

When applied to the description of hadrons, QCD is a highly complex nonperturbative theory that can¬ 
not be solved directly. Numerical simulations of a space-time discretized version of QCD, lattice QCD 
are promising, but still have limitations. Continuum methods using the Bethe-Salpeter/Dyson- 
Schwinger equations leave out some possibly important features of QCD, and so far do not 

describe all decay properties of the meson spectrum and in particular of the heavy-light sector. More¬ 
over, both lattice QCD and Dyson-Schwinger calculations are performed in Euclidean and not in the 
physical Minkowski space. This simplification allows the calculation of mass spectra, but it is not 
clear how the obtained amplitudes can be reliably related to the physical ones. 

This work is inserted in a larger project which proposes to implement the following: 1) a cal¬ 
culation of quark-antiquark bound states in Minkowski space; 2) dynamical quark mass generation, 
due to the self-interaction of quarks calculated consistently from the quark-quark interaction kernel 
and satisfying chiral symmetry constraints; 3) a confining term in the interaction kernel that in the 
nonrelativistic limit reduces to a linear potential. The parameters of the interaction kernel will be 
determined by fits of the meson spectra, and are constrained by the quark mass functions obtained in 
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Figure 1. Diagrammatic representation of the ICSE. The cross indicates that particle 1 is on its mass-shell. P is 
the total momentum and Ejcs the vertex function of the bound-state. 


lattice QCD. This will also provide information on the Lorentz structure of the confining interaction, 
which is not well known at present. The dressed electromagnetic quark current will be calculated 
consistently, using the same kernel. 

Once the quark current is known, many more applications become possible, such as hadron form 
factors and decay properties. We single out the calculation of hadronic contributions to light-by-light 
scattering, which is the most uncertain part of the theoretical prediction of the muon’s anomalous 
magnetic moment. This uncertainty needs to be reduced in order to find out if the observed 3cr 
discrepancy between theory and measurement is a sign of new physics beyond the Standard Model 


2 Pseudoscalar meson spectra with the one-channel CST equation 

We use the Covariant Spectator Theory (CST) 101, which can be viewed as a reorganization of the 
complete Bethe-Salpeter equation (BSE), to develop a dynamical quark model in Minkowski space 
that can describe the structure and the mass spectrum of both heavy and light quark systems. 

The simplest version of the CST equation, called the one-channel CST equation (ICSE) and shown 
diagrammatically in Eig. [T] is obtained when only the positive-energy pole contribution from the 
heavier quark propagator in the energy loop integration is kept. 

An important property of the ICSE is that, unlike the BSE in ladder approximation, it has a smooth 
nonrelativistic limit, and therefore in calculations of heavy quarkonia its results can be compared 
directly to the ones obtained with the Schrodinger equation (SE). 

To find an appropriate way to treat a linear confining interaction in momentum space, and to test 
the required numerical methods, in |01 we used the SE to calculate the two-body bound states of a 
linear potential. Since in S waves these solutions are known exactly, this is an ideal test case. It 
turned out that the singularities arising in the kernel can be treated efficiently, yielding accurate and 
stable numerical solutions. In this work, the same methods were applied for the first time to solve the 
relativistic ICSE. 

The kernel employed in our calculations with the ICSE consists of a covariant generalization 
of the linear (L) confining potential used in |0], a one-gluon exchange (OGE), and a constant (C) 
interaction, 

^{p, k) = ^dp, k) + ^oge(p, k) + ^c(p, k), (1) 
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Figure 2. Mass spectra computed with the ICSE for the bb, cc quarkonia, and for the heavy-light bs, bujbd, be, 
ciilcd, sujsd pseudoscalar states. The horizontal bars represent the experimental 7^^ = 0 states, the centers 
of the crosses and open squares indicate the calculated states for models 1 and 2, respectively. The states fitted 
to the data are those inside the shaded regions, whereas the remaining states inside the dashed rectangles are 
predictions. 


Table 1. Fitted parameters for kernel models 1 and 2. 



cr [GeV^] 

a 

C [GeV] 

y 

Model 1 

0.22 

0.38 

0.337 

8.60 X 10-' 

Model 2 

0.21 

0.37 

0.311 

0.38 X 10-^ 


q - p-kh the transferred four-momentum, and Ep^ = (m^ + is the energy of the heavier quark 
with mass nii. The mixing parameter y allows to dial continuously between a scalar-plus-pseudoscalar 
structure, suggested recently in 111] due to chiral-symmetry constraints, and a vector structure, while 
preserving the same nonrelativistic limit. The precise Lorentz structure of the confining interaction is 
not known, and by fitting the y parameter from the mesonic spectra some further information can be 
gained. The three coupling strengths cr, a and C, were treated as free model parameters. 

For this preliminary work, we used constant masses for the constituent quarks, with the values 

mb = 4.793 GeV, = 1.530 GeV, m, = 0.400 GeV, and = mj = 0.258 GeV . (3) 

The ICSE—with retardation included—requires regularization. We used a Pauli-Villars regulariza¬ 
tion with a cut-off parameter proportional to the average mass of the two constituent quarks, m: 
A = 1.7m. 

The pseudoscalar meson spectra were fitted with two distinct models with different Lorentz 
structure of the confining interaction; whereas model 2 has exactly the equally weighted scalar- 
pseudoscalar component, mixed with a vector structure, as given in Eq. (|2]i, model 1 omits the pseu¬ 
doscalar part. The obtained results are given in Eigure|2]and Table [T] 
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We find that good fits of the pseudoscalar meson spectra can be obtained with either model, and 
that the predictions for the mesons containing lighter quarks are also very reasonable. The observed 
slightly larger difference between the model predictions and the data for the lightest systems is in 
part due to a pole in the propagator of the lighter quark whose importance increases with decreasing 
bound-state mass, but which is omitted in the ICSE. This pole will be included in future work. 

It turns out that both model fits favor essentially pure Lorentz scalar structures of the confining 
interaction. It will be interesting to see if this feature persists once the fits are extended to include also 
vector mesons. 

The inclusion of a pseudoscalar part in the confining interaction leads to a marginally better fit. 
However, because the heavier mesons are almost nonrelativistic systems for which the pseudoscalar 
interaction is strongly suppressed, the difference between models 1 and 2 becomes noticable only for 
the lighter states. In these cases, on the other hand, it may become necessary to replace the ICSE by a 
more complete version of the CST bound-state equations, namely the two-channel (2CSE) or—in the 
case of the pion—the complete four-channel spectator equation (4CSE) i^. 

3 Conclusions and Outlook 

The numerical techniques developed so far allowed us to solve the ICSE and perform the first fit with 
a realistic linear-plus-Coulomb-plus-constant kernel to a wider range of quarkonia and heavy-light 
pseudoscalar states. 

As an immediate continuation of this work, we plan to extend the fit to other mesons (most impor¬ 
tantly vector mesons), and to investigate the effect of replacing the constant constituent quark masses 
by quark mass functions that are calculated in a fully self-consistent fashion from the quark-antiquark 
kernel. The next goal will be to solve the 2CSE and 4CSE to obtain reliable solutions also for the light 
sector. 

Meanwhile, based on the results determined so far, we can already state that the calculational 
methods we developed to treat a linear confining interaction in momentum space and tested in the 
nonrelativistic SE, also work in the relativistic case of the ICSE in Minkowski space. These tech¬ 
niques will be useful in future studies of the mesonic bound-state problem using the CST. 
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